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NOTATION 

dp/dn  =  pressure  gradient 
f  =  Coriolis  parameter 
K  =  coefficient 

^  =  angle  between  storm  track  and  coastline 

p  =  retardation  factor 
p^  =  standard  pressure 
P(n,''^)  =    probability  that  n  storms  occur  in  t  years 
P(V  ^  vin)  =  probability  that  the  wind  speed  V  in  n  storms  is  less  than 
P(V      V I ''^)  =  probability  that  the  wind  speed  V  in      years  is  less  than  \ 

r  =  distance  from  storm  center 

R  =  radius  of  maximum  wind  speeds 

s  =  hurricane  translational  speed 

t  =  time,  in  hrs. 

V      =  maximum  gradient  windspeed 
gx  ^  ^ 

V(z,R)  =  maximum  velocity  at  z  meters  above  ground  at  the  radius  of 
maximum  winds 

V(z,r)  =  maximum  velocity  at  z  meters  above  ground  at  the  radius  r 
from  the  eye 

V(z,r,®)  =  velocity  at  z  meters  above  ground  at  radius,  r  and  angle 
^  from  OM 

z  =  height  above  ground  or  water  surface,  in  meters 
"  =  coefficient 

max  ~  pressure  drop  from  periphery  of  hurricane  to  its  center 
p 

=  air  density 

0 

=  angle  between  lines  ON  and  OM 

(J) 

=  angle  between  storm  track  and  coastline 
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ABSTRACT 


A  Monte  Carlo  simulation  technique  Is  used  to  obtain  estimates  of 
hurricane  wind  speeds  along  the  Gulf  and  East  Coasts  of  the  United 
States.    The  paper  describes  the  sources  of  data,  the  probabilistic 
models  for  cllmatologlcal  characteristics  of  hurricanes,  and  the  physical 
models  for  the  hurricane  wind  speed  field  used  In  the  estimations. 
Estimated  values  of  fastest-mile  hurricane  wind  speeds  at  10  m  above 
ground  In  open  terrain  at  the  coastline  and  at  200  km  Inland  are  given 
for  various  mean  recurrence  Intervals.    The  estimated  hurricane  wind 
speeds  were  found  to  be  best  fitted  by  Welbull  distributions  with  tall 
length  parameters  y      U.    Estimates  are  given  of  various  errors  Inherent 
In  the  estimated  values  of  the  hurricane  wind  speeds.    Owing  to  uncertain 
ties  with  respect  to  the  applicability  of  the  physical  models  used 
In  this  work  to  locations  north  of  Cape  Hatteras,  estimated  hurricane 
wind  speeds  given  for  these  locations  should  be  viewed  with  caution. 

Key  Words:     Buildings  (codes);  climatology;  hurricanes;  statistical 

analysis;  structural  engineering;  tropical  cyclones;  wind 
(meteorology) . 


Facing  page;    Damage  caused  by  Hurricane  Camille, 
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1 .  INTRODUCTION 


It  has  previously  been  shown  that  predictions  of  extreme  wind  speeds 
in  hurricane-prone  regions  cannot  in  general  be  based  upon  the  statis- 
tical analysis  of  the  largest  annual  wind  speeds  recorded  at  a  given 
site  [10,  p.  84].     For  this  reason,  estimates  of  extreme  hurricane 
wind  speeds  at  a  site  are  commonly  obtained  by  indirect  methods  [1, 
4,  11,  12]   from  (1)  statistical  information  on  the  climatological 
characteristics  of  hurricanes,  and  (2)  a  physical  model  of  the  hurri- 
cane wind  structure.     The  climatological  characteristics  of  hurricanes 
include:     (a)  rate  of  hurricane  occurrence  in  any  given  region;  (b) 
difference  between  atmospheric  pressures  at  the  center  and  at  the 
periphery  of  the  storm;  (c)  radius  of  maximum  wind  speeds;   (d)  speed 
of  storm  translation;  (e)  direction  of  storm  motion,  and  (f)  crossing 
point  coordinate  along  the  coast  or  on  a  line  normal  to  the  coast. 
The  physical  model  of  the  hurricane  wind  includes  assumptions  on: 

(a)  dependence  of  surface  wind  speeds  upon  difference  between  atmos- 
pheric pressures  at  center  and  periphery  of  the  storm,  radius  of 
maximum  wind  speeds,   speed  of  translation,   latitude,  and  position 

of  the  point  being  considered  with  respect  to  center  of  the  storm; 

(b)  storm  decay  as  the  storm  travels  over  land  and  its  supply  of 
energy  in  the  form  of  warm  moist  air  from  the  ocean  surface  is 
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thus  cut  off;  (c)  reduction 
and  (d)  ratios  between  wind 


of  wind  speeds  due  to  friction  over  land; 
speeds  averaged  over  various  time  intervals. 


Given  the  available  statistical  information  on  the  climatological  charac- 
teristics of  hurricanes,  probabilistic  models  are  estimated  for  each  of 
these  characteristics.     A  large  number  of  hurricanes  is  then  generated 
by  Monte  Carlo  simulation  (random  sampling)  from  these  probabilistic 
models,  as  will  be  explained  subsequently  in  the  report.     The  largest 
wind  speeds  occurring  within  each  of  these  hurricanes  at  the  site  of 
concern  are  the  data  from  which  the  cumulative  distribution  function 
of  the  hurricane  wind  speeds  at  that  site  is  then  estimated.  This 
procedure  was  first  developed  by  Russell  [9],  who  also  developed  the 
computer  program  used  in  this  work. 

The  purpose  of  this  paper  is  to  present  estimates  of  hurricane  wind 
speeds  on  the  Gulf  and  East  Coasts  of  the  United  States  based  on  the 
procedure  just  outlined  and  on  data  taken  from  References  3  and  5.  The 
probabilistic  and  physical  models  used  to  obtain  these  estimates  will 
be  described,  and  an  analysis  will  be  made  of  errors  inherent  in  the 
estimates.     Estimates  of  hurricane  wind  speeds  in  California  [13]  are 
not  treated  in  this  report. 


Facing  page:    Damage  to  solar  collectors,  caused  by  Hurricane  Frederick. 
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2.     PROBABILISTIC  AND  PHYSICAL  MODELS 
2 . 1     PROBABILISTIC  MODELS 

The  following  probabilistic  models  were  used: 

1.  The  hurricane  occurrence  is  described  by  a  constant  rate 
Poisson  process. 

2.  The  probability  distribution  of  the  pressure  difference  between 
center  and  periphery  of  storm,  APjjj^jj,  is  lognormal.     To  eliminate 
values  of  APjjj^j^  judged,  in  light  of  historical  data  [5],  to  be 
unrealistically  high,  the  distribution  is  censored  so  that 
AP„o^     <  101.6  mm  (4.00  in)  of  mercury.     (Note  that  APm-,^ 

LUd^  lUdA 

101.6  mm  corresponds  to  the  lowest  atmospheric  pressure  ever 
recorded  worldwide  [14].)    The  effect  of  this  censoring  will  be 
commented  upon  subsequently  in  this  report. 

3.  The  probability  distribution  of  the  radius  of  maximum  wind  speeds, 
R,  is  lognormal.  This  distribution  is  censored  so  that  8  km  <  R  < 
100  km  to  avoid  unrealistically  "tight"  or  "broad"  storms  [5]. 
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4.  The  average  correlation  coefficient  between  R  and  A Pmax 
approximately  -0.3  (see  Reference  5,  pp.  68  and  69).     All  other 
climatological  characteristics  are  assumed  to  be  statistically 
independent . 

5.  The  probability  distribution  of  the  speed  of  translation,  s,  is 
normal.  This  distribution  is  censored  so  that  2  km/hr  <  s  <  65 
km/hr  [5] . 

6.  The  cumulative  distribution  function  of  the  hurricane  crossing 
point  along  the  coast  is  a  curve  matching  the  historical  data 
as  recorded  in  Reference  5.  Separate  distribution  curves  are 
defined  for  entering,  exiting,  upcoast  heading,  and  downcoast 
heading  storms.     The  length  of  coast  being  considered  for 
entering  and  exiting  storms  includes  a  470  km  segment  downcoast 
of  the  location  under  investigation,  and  a  370  km  segment  up- 
coast of  that  location,   the  influence  of  storms  not  crossing 
these  segments  being  negligible  for  practical  purposes. 

The  downwind  and  upwind  segments  are  unequal  because,  owing 
primarily  to  the  presence  of  the  translational  speed  of  the 
storm,  the  hurricane  wind  structure  is  asymmetrical.  As  an 
example,  the  distribution  curve  assumed  for  entering  storms 
affecting  coastal  milepost  400  n.  mi.  (740  km)  [see  Figure  1] 
is  shown  in  Figure  2. 

7.  The  crossing  points  on  a  line  normal  to  the  coast  are  assumed 
to  be  uniformly  distributed  between  110  km  and  280  km  from  the 
coast  in  the  case  of  upcoast  heading  storms,   and  between  110  km 
and  360  km  in  the  case  of  downcoast  heading  storms. 

8.  The  cumulative  distribution  functions  of  heading  (direction  of 
storm  translation)  for  entering  storms  are  curves  matching  the 
historical  data  recorded  in  Reference  5.    As  an  example.  Figure  3 
shows  the  distribution  of  the  storm  headings  about  their  mean 
value  for  milepost  400  n.  mi.  (740  km).     For  exiting  storms  and 
for  storms  crossing  a  line  normal  to  the  coast,  the  distributions 
are  assumed  to  be  uniform  between  +  30°  of  the  mean.     The  mean 
values  of  the  storm  headings  are  obtained  from  the  data  of 
References  3  and  5.     In  all  cases  it  was  assumed  that  the  path 

of  the  storm  motion  is  a  straight  line. 

2.2     PHYSICAL  MODELS 

2.2.1     Maximum  Gradient  Wind  Speed 

The  maximum  gradient  wind  speed  can  be  written  as 

V         _  Rf  +         dp^l/2  (1) 
Vgx  .      -J  ^  ^- 
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Figure  1.    Locator  map  with  coastal  distance  intervals  marked 
in  nautical  miles  (1  n.  mi.   =  1.9  km)  [5]. 
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(see  Reference  10,  p.  18),  where  f  =  Coriolis  parameter,  dp/dn  =  pressure 
gradient,   P  =  air  density,  and  R  =  radius  of  maximum  wind  speeds.     It  is 
reasonable  to  assume 

^  =  «^Pmax  (2) 
an 

where  «     is  a  coefficient  determined  empirically.     It  then  follows  from 
Equations  1  and  2 

Vgx  =  K    V  Apmax  -  ^  (3) 

1/2 

where  the  notation  K  =  [R  "/P]        is  used.     Values  of  K  consistent  with 
those  determined  empirically  by  the  Hydrometeorological  Branch  (National 
Weather  Service,  National  Oceanic  and  Atmospheric  Administration)  and 
used  in, Reference  8  vary  approximately  between  6.97  m/s/mm   ,  (77.5 
mph/in^^  )  at  latitude  23°N  to  6.93  m/s/mm^'     (77.0  mph/in-""^  )  at 
latitude  45°  N.     [These  values  may  be  used  in  conjunction  with  the 
assumption  that  the  pressure  at  the  periphery  of  the  storm  is  p^  =  765 
mm  (29.77  in).] 

2.2.2    Wind  Speeds  at  10  m  Above  the  Ocean  Surface 

The  maximum  wind  speed  at  10  m  above  the  ocean  surface,  averaged  over 
10-min,   is  assumed  to  be  given  by  the  empirical  relation  [8] 


V(z=10,R)  =  0.865  V      +  0.5  s  (4) 

gx 

Let  the  center  of  the  storm  be  denoted  by  0,  and  let  a  line  OM  be 
defined  that  makes  an  angle  of  115°  clockwise  from  the  direction  of 
motion  of  the  storm.     The  10-min  wind  speed  at  10  m  above  the  ocean 
surface  at  a  distance  r  from  0  along  the  line  OM  is  denoted  by  V(z=10, 
r).     The  ratio  V(z=10,r)/V(z=10,R) ,  is  assumed  to  depend  on  r  as  shown 
in  Figure  4  [8] .     Let  now  the  angle  between  a  line  ON  and  line  OM  be 
denoted  by  6.     The  10-min  wind  speed  V(z=lO,r,0)  at  10m  above  the 
surface  at  a  distance  r  from  the  storm  center  on  line  ON  is  assumed 
to  be  given  by  the  expression  [8] 

V(z=10,r,e)  =  V(z=10,r)  -^(1  -  cos  6)  (5) 

The  wind  velocity  vector  has  a  component  directed  toward  the  center  of 
the  storm.     The  angle  between  that  vector  and  the  tangent  to  the  circle 
centered  at  0  is  assumed  to  be  between  0°  and  10°  in  the  region 
0  <  r  <  R,  between  10°  to  25°  in  the  region  R  1  r  <  1.2R,  and  25°  in 
the  region  r  1  1 .2R  [8] . 


8 


2.2.3  Storm  Decay 


The  storm  decay  was  assumed  to  result  from  a  decrease  with  time  of  the 
difference  between  pressure  at  the  center  and  at  the  periphery  of  the 
storm,  given  by  the  following  relation 

^p(t)  =  -  0.02  [1+sin  f]t  (6) 

where  t  =  travel  time  in  hours,  ^p(t)  and  ^p^^^      are  given  in  inches, 
and    't'  =  angle  between  coast  and  storm  track        <  <t>  <  180").    As  shown 
subsequently  in  this  paper,  results  based  on  this  model  were  compared 
with  results  based  on  Malkin's  model,  which  is  consistent  with  measure- 
ments reported  in  Reference  7.     This  comparison  appears  to  support  the 
validity  of  Equation  6  as  a  conservative  approximation. 

2.2.4  Reduction  of  Wind  Speeds  Due  to  Friction  Over  Land 

The  factor  for  reducing  overwater  to  overland  surface  wind  speeds  used 
in  Refej^ence  8  is  0.78  for  wind  speeds  over  37  m/s.     In  Reference  2,  the 
ratio  V  (10)/v"(10),  where  V(10)  =  10-min  mean  speed  at  10  m  above  the 
surface,  and  the  superscripts  w  and  ^  indicate  "overwater"  and  "overland' 
respectively,  is  given  as 

V^(IO)  ^  1  (7) 

10 

V^(IO)  P 

where  the  retardation  factor  p  =  0.83  and  the  roughness  length  Zq  = 
0.005  m,  so  that  V  (10)/V^(10)  =  0.79,  i.e.,  this  ratio  has  almost 
exactly  the  same  value  as  in  Reference  8.     Note  that  Equation  7  is  based 
on  similarity  relations  applicable  to  storms  in  which  the  gradient  wind 
may  be  regarded  as  geostrophic,  so  that  its  applicability  to  hurricane 
winds  is  uncertain,  at  least  in  the  region  of  maximum  wigd  speeds.  In 
this  work  it  was  therefore  assumed  conservatively  that  V  (10)/V^(10) 
=0.85. 

2.2.5  Dependence  of  Wind  Speeds  Upon  Averaging  Time 

Research  results  concerning  the  ratios  of  hurricane  wind  speeds  corres- 
ponding to  different  averaging  times  do  not  appear  to  be  available. 
In  the  absence  of  such  results,  it  was  judged  reasonable  to  use  wind 
speed  ratios  obtained  for  non-tropical  storms  by  Durst  which  are  summa- 
rized, e.g.,  in  Reference  10,  p.  62.     For  example,  the  ratio  between 
the  fastest  one-minute  and  the  fastest  10-minute  speed  at  10  m  above 
open  terrain  was  assumed  to  be  1.18. 


Facing  page:     Trees  uprooted  by  typhoon  in  Manila  (The  Philippines), 
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3.     PROBABILITIES  OF  OCCURRENCE  OF  HURRICANE  WIND  SPEEDS 


3.1     THEORETICAL  APPROACH 

To  estimate  probabilities  of  occurrence  of  hurricane  wind  speeds,  one 
thousand  hurricanes  were  assumed  to  hit  the  area  adjoining  each  of  the 
56  sites  being  investigated.     (Recall  that  this  area  was  defined  pre- 
viously in  the  section  "Probabilistic  Models",  items  6  and  7.)  Four 
categories  of  storms  were  considered:  entering  hurricanes  (i.e.,  hurri- 
canes moving  inland),  exiting  hurricanes  (i.e.,  hurricanes  headed 
toward  the  ocean  after  moving  over  land),  upcoast  heading  hurricanes 
(i.e.,  hurricanes  whose  center  remains  offshore  and  moving  upcoast), 
and  downcoast  heading  hurricanes.     For  each  location,  the  ratios  between 
the  number  of  hurricanes  belonging  to  each  of  these  categories  and  the 
total  number  of  hurricanes  affecting  that  location  were  estimated  from 
the  information  given  in  References  3  and  5. 

The  climatological  data  for  each  of  the  one  thousand  hurricanes  were 
determined  from  the  respective  probability  distributions  by  Monte  Carlo 
simulation  (random  sampling).     Associated  with  each  hurricane  is  a  wind 
field,  and  a  set  of  wind  speeds  at  the  site  of  interest  which  depend 
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upon  the  position  of  the  moving  hurricane  center  with  respect  to  the 
site.     These  wind  speeds  were  calculated  for  a  sufficiently  large 
number  of  such  positions  (6  or  7,  depending  on  track).     The  largest 
among  these  speeds  represents  the  maximum  wind  speed  caused  by  the 
hurricane  at  the  site.     A  set  of  one  thousand  such  wind  speeds  is 
thus  obtained,  which  is  used  as  the  basic  set  of  data  for  the  esti- 
mation of  the  probability  of  occurrence  of  hurricane  wind  speeds 
at  the  site.     (This  procedure  was  applied  to  estimate  extreme  wind 
speeds  without  regard  to  direction,  as  well  as  extreme  wind  speeds 
blowing  in  specified  directions.)     The  wind  speeds  are  ranked  by  magni- 
tude, and  the  probability  of  occurrence  of  the  i-th  ranking  wind 
speed,  v^,  in  a  set  of  m  wind  speeds  (m  =  1,000  in  the  work  reported 
herein)  is  obtained  as  follows. 

Let  the  probability  that  the  wind  speed  in  any  one  storm  is  less  than 
V  be  denoted  by  F^.  The  probability  that  the  highest  wind  speed  V  in 
n  storms  is  less  than  v  can  be  written  as 

n  ^  (8) 


P(V  <  v|n)     =  (F^)"  =  F 


V 


Let  the  probability  that  V  <  v  in  t  years  be  denoted  by  P(V  <  v,  x). 
It  is  possible  to  write 


P(V  <  V,  t)     =    Z  p(v  <  v|n)  P(n,  t)  (9) 
n=0 


where  P(n,T)  denotes  the  probability  that  n  storms  will  occur  in  t 
years.     If  a  Poisson  process  is  assumed  to  describe  P(n,T),  and  if  use 
is  made  of  Equation  9, 

°°  (  VT  )  ^      — VT 

P(V  <  V,  t)     =     Z     F'^  ® 
n=0    ^  n! 

"     (vtF  )^ 
=    e"  ^"^  Z  ^ 

'^=0  n! 

_    VT  VT^ 

=    e  e  V 


VT 


=  e 


(l-Fv^  (10) 


where  the  Poisson  parameter  ^'  =  annual  rate  of  occurrence  of  hurricanes 
in  the  area  of  interest  for  the  site  being  considered.  For  t  =  1, 
P(V  <  V,  t)  =  probability  of  occurrence  of  wind  speeds  less  than 
V  in  any  one  year . 
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Consider  now  the  wind  speed  .  Its  probability  of  occurrence  can  be 
written  as 


Thus 

P(V  <  VI, 1)  =  e  ^  (12) 

The  estimated  mean  recurrence  interval  of  the  speed  v^^  is  then 

N  i   (13) 

1-e 

3.2     NUMERICAL  RESULTS 
3.2.1    Extreme  Wind  Speeds 

Hurricane  fastest-mile  wind  speeds  at  10m  above  ground  over  open  terrain 
were  estimated  for  58  mileposts  (from  150  through  3000  in  increments  of 
50  -  see  Figure  1)  at  the  coastline,  and  at  25  km,  50  km,  100  km  and 
200  km  inland.     The  results  of  the  calculations  were  smoothed  as  follows. 
For  each  milepost  the  average  was  calculated  of  the  fastest-mile  wind 
speeds  at  that  milepost  and  at  its  two  neighboring  mileposts  [located 
each  at  50  n.mi.  (93  km)  from  it].     This  averaging  was  done,  for  each 
milepost,  at  the  coastline,  as  well  as  at  25  km,  50  km,  100  km,  and 
200  km  inland.     A  least-squares  straight  line  was  fitted  at  each  mile- 
post to  the  five  averages  thus  obtained.     The  ordinates  of  this  line  at 
the  coastline  and  at  200  km  inland  represent  the  values  given  by  the 
curves  of  Figures  5  and  6,  respectively.    Extreme  fastest-mile  speeds 
at  distances  inland  of  less  than  200  km  from  the  coastline  may  be 
obtained  by  linear  interpolation  between  Figures  5  and  6.     It  is  noted 
that  the  values  of  Figures  5  and  6  differ  by  at  most  a  few  percent  from 
the  corresponding  unsmoothed  values  originally  calculated  for  the  coast- 
line locations. 

The  results  shown  in  Figures  5  and  6  represent  smoothed  estimates  of 
hurricane  wind  speeds  without  regard  to  direction.     As  previously 
indicated,  similar  results  were  obtained  for  winds  blowing  in  specified 
directions.     Estimates  of  coastline  fastest-mile  hurricane  wind  speeds 
at  10  m  above  ground  over  open  terrain,  smoothed  by  averaging  results 
obtained  for  three  consecutive  milestone  positions,  are  plotted  in 
Appendix  A  as  a  function  of  direction.     Note  that  the  number  of  storms 
associated  with  strong  winds  blowing  in  specified  directions  is  in 
certain  cases  very  small.     It  was  therefore  not  possible  in  those 
cases  to  estimate  wind  speeds  corresponding  to  relatively  short  mean 
recurrence  intervals.     This  is  reflected  by  the  blank  portions  of 
some  of  the  curves  of  Appendix  A. 
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Errors  in  the  estimation  of  extreme  hurricane  wind  speeds  are  dealt 
with  in  the  subsequent  section.     It  is  noted  here  that  these  errors 
increase  as  the  mean  recurrence  interval  increases.    Therefore,  esti- 
mates of  wind  speeds  corresponding  to  mean  recurrence  intervals 
exceeding  100-yr,  say,  may  contain  significant  errors.  Nevertheless, 
such  estimates  obtained  by  the  procedure  used  in  this  paper  are  of 
interest  as  they  provide  the  only  available  basis  for  assessing  the 
validity  of  wind  load  factors  currently  being  used  for  structural 
design  purposes  in  hurricane-prone  regions. 

3.2.2    Probability  Distributions  of  Hurricane  Wind  Speeds 

The  results  obtained  from  the  calculations  were  used  to  estimate  proba- 
bility distributions  of  hurricane  wind  speeds  (blowing  from  any  di- 
rection).    It  was  found,  by  using  the  probability  plot  correlation 
coefficient  method,  that  in  all  cases  the  best  fitting  distribution 
for  the  1-minute  speeds  was  the  Weibull  distribution  with  tail  length 
parameter  Y  >  4  (for  practical  purposes  Y  =  4  can  be  assumed,  the 
differences  between  results  based  on  Y  =  4  and    Y  >  4  being  negligible). 
As  an  illustration,  the  results  obtained  at  milepost  450  [1-minute 
wind  speeds  at  coastline  over  water,  in  knots  (1  knot  ~  1.9  km/hr)] 
are  plotted  on  Weibull(Ys4)  and  Extreme  Value  Type  I  probability  paper 
on  Figures  7  and  8,  respectively.     It  is  seen  that  the  results  fit 
the  Extreme  Value  Type  I  distribution  poorly,  i.e.,  the  results  do 
not  fit  a  straight  line  on  Fig.  8,  as  they  do  on  Fig.  7. 

3.3    MIXED  PROBABILITY  DISTRIBUTIONS  OF  HURRICANE  AND  NON-HURRICANE 
WIND  SPEEDS 

Mixed  probability  distributions  of  hurricane  and  non-hurricane  fastest- 
mile  (as  opposed  to  1-minute)  wind  speeds  were  estimated  by  using  the 
following  expression: 

P(V  <  V)  =  Pjj  (Vjj  <  v)-Pnh(Vnh  <  V)  (14) 

in  which  Pjj  (Vjj  <  v)  =  cumulative  distribution  of  hurricane  wind  speeds 
Vjj,  and  Pj^jj  ^^NH  ^  ^)  ~  cumulative  distrtibution  of  non-hurricane  wind 
speeds,  Vjj^j.     The  distribution  Pj^  was  determined  as  shown  in  this  report. 
For  the  distribution  Pfjjj,  an  Extreme  Value  Type  I  distribution  was 
assumed.     The  parameters  of  Pjjjj  were  estimated  at  various  locations  from 
maximum  annual  wind  speed  data  that  did  not  include  hurricane  wind 
speeds.     The  results  of  the  calculations  showed  that  the  effect  of  the 
non-hurricane  winds  is  negligible  for  mean  recurrence  intervals  of  the 
order  of  50  years  or  more  (see  Fig.  9).     For  mean  recurrence  intervals 
of  about  20  years,  the  estimated  wind  speeds  that  include  the  effect 
of  non-hurricane  winds  exceed  the  estimated  hurricane  speed  by  about 
5%.    Note  that  these  conclusions  are  not  applicable  north  of  Cape 
Hatteras,  where  non-hurricane  wind  may  control  the  design  at  certain 
locations . 


16 


o 


u 


=3 


®    <S>  <S  <S>      «6     ®      Lf)  (X) 

s>  <s  <s  in     ru  «-i 
in  ru 


ru 


o 
.J 
a. 

> 


CO 

<i 

QQ 
O 

q: 

Q. 


3 


> 


Ui 
O 


u. 


(0 


01 


9 


LA 


9 


ID 


at 


u 

O 
Ui 
K 


17 


o 

bJ 
OL 

z 

Z) 


9 


ru 


9 


ru  "-I 


9 

Lft 


ru 


it) 


nj 


o 
a 


o 


a. 


X 
UJ 


9 


00 


<J3 


18 


o 

t-t 

a. 
ill 

a. 

z 
Q: 

Z3 
t- 
llJ 

a: 


9 

nj 


o 
a. 

> 


« 

M 

O 

0. 


a. 
> 


t- 
X 


9 

9 

• 

• 

• 

• 

• 

• 

* 

1^ 

U) 

(U 

1 

• 

a- 

B 

• 

W 

w 

73 

<U 

(U 

Cu 

0) 

• 

-a 

wi 

(U 

iH 

•H 

• 

(3 

1 

4-* 
(0 

0) 

4J 

(0 

1  ft 

Ul 

IM 

O 

CO 

C 

o 

•H 

■U 

3 

^ 
•H 

(0 

•H 

• 

t3 

t'^ 

LT) 

■U 

•H 

iH 

•H 

.o 

ns 

• 

o 

in 

M 
O. 

H 

0) 

• 

a. 

L/) 

H 

(U 

iH 

> 

• 

a 

■U 

X 

00 
•H 


19 


Facing  page:    Damage  caused  by  Cyclone  Taray  in  Darwin,  Australia. 
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4.     ESTIMATION  ERRORS 


The  errors  in  the  estimation  of  hurricane  wind  speeds  by  the  procedure 
used  in  this  work  may  be  divided  into  four  categories:     (1)  sampling 
errors ,  due  to  (a)  the  limited  size  of  the  data  sample  used  in  making 
statistical  inferences  on  the  climatologogical  characteristics  of 
hurricanes  (in  the  United  States  these  data  samples  correspond  to 
lengths  of  record  of  75  to  100  years),  and  (b)  the  limited  number  of 
hurricanes  generated  by  the  Monte  Carlo  simulation;  (2)  probabilistic 
modeling  errors,  due  to  the  imperfect  choice  of  the  distribution 
functions  to  which  the  climatological  data  are  fitted  (e.g.,  assuming 
that  a  lognormal  distribution  holds  when  in  fact  a  normal  distribution 
would  be  more  appropriate);  (3)  observation  errors,  due  to  the  imper- 
fect measurement  or  recording  of  the  true  values  of  the  climatological 
characteristics;  and  (4)  physical  modeling  errors,  due  to  the  imperfect 
representation  of  the  dependence  of  the  wind  speed  upon  the  various 
climatological  characteristics  and  micrometeorological  parameters. 

4.1     SAMPLING  ERRORS 

As  shown  in  some  detail  in  Reference  1,  the  standard  deviation  of  the 
errors  associated  with  the  limited  size  of  the  climatological  data 
samples  (corresponding  to  about  75  to  100  years  of  record)  is  of  the 
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order  of  6%  to  10%  of  the  estimated  wind  speeds,  while  the  standard 
deviation  of  the  errors  due  to  the  limited  number,  m,  of  hurricanes 
used  in  the  simulation  is  of  the  order  of  a  few  percent  for  m  j,  1,000. 
In  some  cases,  say,  at  Jacksonville,  Florida,  the  poorly  understood 
and  seemingly  anomalous  occurrence  histories  may  lead  to  larger  errors 
than  are  typical  of  the  model  overall.     Part  of  the  sampling  errors  due 
to  the  limited  number  of  hurricanes  used  in  the  simulation  is  elimi- 
nated by  the  use  of  the  smoothing  procedure  described  previously. 


4.2     PROBABILISTIC  MODELING  ERRORS 

To  assess  the  effect  of  various  probabilistic  models  used  in  the  simu- 
lation, calculations  were  carried  out  at  mileposts  400,  1450  and  2100 
(see  Figure  1)  under  the  following  assumptions: 

a.  Distribution  of  aP        is  normal.     This  assumption  results  in  a 

max 

decrease  of  the  50-yr,  100-yr,  and  400-yr  wind  speeds  with  res- 
pect to  the  corresponding  values  based  on  the  assumption  of 
lognormality  by  about  5%,  5%,  and  8%,  respectively. 

b.  Distribution  of  R  is  normal.     This  assumption  results  in  an 
increase  of  the  50-yr,  100-yr,  and  400-yr  wind  speeds  with 
respect  to  their  corresponding  values  based  on  the  assumption 
of  lognormality  by  about  1%,  2%,  and  6%,  respectively. 

c.  Distributions  of  APmax         ^  both  normal.     This  assumption 
results  in  a  decrease  of  the  50-yr,  100-yr,  and  400-yr  wind 
speeds  with  respect  to  their  corresponding  values  based  on  the 
assumption  that  APmax  ^^'^  ^          both  lognormally  distributed 
by  about  3%,  3%,  and  4%,  respectively. 

d.  Distribution  of  s  is  lognormal.     The  differences  between  esti- 
mates of  wind  speeds  based  on  this  assumption  on  the  one  hand, 
and  on  the  assumption  that  s  is  normally  distributed  on  the 
other  hand,  are  negligible  for  practical  purposes. 

e.  Distribution  of     ap        is  not  censored.     The  censoring  of  the 

^  max  ° 
distriibution  of  APmax  ^ ^^max      101 -6  mm)  reflects  the  writers' 
belief  that  this  upper  bound  is  dictated  by  physical  consider- 
ations, and  that  an  unbounded  distribution  would  thus  be 
physically  unrealistic.     Nevertheless,  calculations  were  carried 
out  for  five  stations  in  which  it  was  assumed  that  the  lognormal 

distribution  of  ap        is  unbounded.     It  was  found  that  the 

max 

estimated  hurricane  wind  speeds  and  their  cumulative  distribution 
functions  were  for  practical  purposes  the  same  regardless  of 
whether  the  condition  .p  101.6  mm  was  assumed  or  not. 
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A  final  observation  with  regard  to  the  choice  of  probabilistic  models 
concerns  the  representation  of  the  hurricane  occurrences  as  a  Poisson 
process.  It  can  be  shown  that  this  model  could  be  replaced,  e.g.,  by 
a  geometric  probability  model,  without  any  significant  effect  upon  the 
estimated  results,  particularly  for  winds  with  mean  recurrence  inter- 
vals of  the  order  of  20  years  or  more. 

4.3  OBSERVATION  ERRORS 

To  assess  the  effect  of  possible  errors  in  the  measurement  (or  re- 
cording) of  the  parameters  ^P„„„  and  R,  calculations  were  carried  out 
at  mileposts  400,  1450,  and  2100  (Figure  1)  in  which  it  was  assumed 
that  the  mean  and/or  standard  deviation  of  ^P^iax  ^  larger 

than  those  estimated  from  the  data  of  Reference  5.     Increasing  both 

the  mean  and  the  standard  deviation  of  ^p        by  a  factor  of  1.1  resulted 

max 

in  an  increase  in  the  estimated  values  of  the  50-yr,  100-yr,  and  400-yr 
winds  of  about  5%  at  milepost  400,  8%  at  milepost  1450,  and  10% 
at  milepost  2100.     Leaving  the  mean  unchanged  and  increasing  the  stan- 
dard deviation  by  a  factor  of  1.1  resulted  in  an  increase  of  the  esti- 
mated speeds  of  the  order  of  a  few  percent.     The  effect  of  increasing 
or  decreasing  the  mean  and/or  the  standard  deviation  of  the  radii  of 
maximum  wind  speeds,  R,  by  a  factor  as  high  as  1.5  was  of  the  order 
of  a  few  percent  at  most< 

As  previously  indicated,   the  data  used  in  this  report  were  obtained 
from  Ref .  5.  The  estimation  of  hurricane  wind  speeds  at  various  loca- 
tions could  presumably  be  improved  by  using  additional  sources  of 
climatological  information,   e.g.,   those  available  in  the  archives 
of  the  Corps  of  Engineers.     The  writers  believe  that  this  is  the  case 
particularly  for  South  Texas,  Louisiana,  southwest  Florida  and  Jackson- 
ville . 

Note  that  estimated  hurricane  wind  speeds  are  lower  at  mileposts  1000 
and  1750  (Figures  1,  5,  and  6)  than  in  the  regions  adjacent  to  these 
mileposts.     This  may  be  due  to  the  relation  between  configuration  of 
the  coast  and  the  hurricane  paths.     However,   it  may  be  that  hurricanes 
at  and  near  these  mileposts  were  underreported  in  the  past.     It  appears 
therefore  prudent  to  increase  the  wind  speeds  calculated  at  mileposts 
1000  and  1750  by  about  5  mph  or  10  mph.     Finally,   it  was  found  that 
the  storm  decay  overland  was  very  weak  at  mileposts  650  and  2250. 
This  is  due,  at  least  in  part,  to  the  configuration  of  the  coast  at 
these  mileposts,  although  the  spurious  occurrence  of  strong  simulated 
storms  may  also  have  contributed  to  this  result. 

4.4  PHYSICAL  MODELING  ERRORS 

The  coefficient  K  in  Equation  3,   the  coefficients  of  Equation  4,  and 
the  dependence  of  wind  speeds  upon  radius  represented  in  Figure  4, 
have  been  determined  empirically  by  the  National  Weather  Service  on 
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the  basis  of  careful  correlations  of  pressure  and  wind  speed  measure- 
ments.    Clearly,  the  corresponding  physical  model  of  the  hurricane  wind 
field  is,  nevertheless,  imperfect.     It  is  the  writers'  belief,  based, 
e.g.,  on  information  from  Reference  6,  that  the  standard  deviation  of 
the  errors  inherent  in  this  model  is  of  the  order  of  5%  to  10%. 

It  was  noted  previously  that  the  coefficient  K  assumed  in  this  work  may 
be  used  in  conjunction  with  the  value  p^  =  756  mm  (29.77  in).     To  check 
the  effect  on  the  calculated  results  of  the  assumption  p^  =  760  mm 
(29.92  in),  calculations  based  on  this  assumption  were  carried  out  at 
mileposts  400,  1450,  and  2100.     It  was  found  that  differences  between 
results  based  on  p^  =  756  mm  and  p^  =  760  mm  were  of  the  order  of  3%. 

It  was  also  noted  that  results  based  on  the  model  for  storm  decay  used 
in  this  work  (Equation  6)  were  compared  with  results  based  on  Malkin's 
model,  which  is  consistent  with  measurements  reported  in  Reference  7. 
The  differences  between  the  two  sets  of  results  were  in  all  cases  small 
(of  the  order  of  +  2%). 

In  estimating  wind  speeds  over  water  immediately  off  the  shoreline  it 
was  assumed  that  the  roughness  length  parameter  corresponds  to  "flow 
over  water"  conditions  in  all  directions  (this  will  be  referred  to  as 
"surface  friction  hypothesis  I").     In  reality,  each  location  immediately 
offshore  is  affected  not  only  by  hurricane  winds  which  blow  from  the 
ocean  (and  for  which  "flow  over  water"  conditions  prevail),  but  also 
by  hurricane  winds  blowing  from  inland,  for  which  the  assumption  that 
a  roughness  length  corresponding  to  "flow  over  open  terrain"  conditions 
may  have  to  be  used.     The  hypothesis  that  the  roughness  length  depends 
upon  whether  the  winds  flow  from  the  ocean  or  from  land  will  be  re- 
ferred to  as  "surface  friction  hypothesis  II".     Whether  or  not  "surface 
friction  hypothesis  II"  is  warranted  is  by  no  means  always  certain: 
indeed,  near  the  shoreline  the  terrain  is  likely  to  be  flooded  during 
the  occurrence  of  hurricanes.     Nevertheless,  it  is  of  interest  to  check 
the  differences  between  results  based  on  surface  friction  hypotheses  I 
and  II.     Calculations  carried  out  at  mileposts  400,  850,  1450,  2100,  and 
2800  showed  that  estimated  wind  speeds  corresponding  to  hypothesis  II 
were  lower  than  those  based  on  hypothesis  I  by  an  amount  generally  not 
exceeding  3%. 

Local  wind  intensification  due  to  the  presence  of  rainbands  occur  at 
all  locations  south  of  the  30°  north  parallel.     These  intensifications 
might  have  to  be  be  accounted  for  by  increasing  the  wind  speeds  cal- 
culated herein  by  as  much  as  10  mph  (4.47  m/s)  or  so. 

As  previously  noted,  the  physical  models  used  in  this  report,  and 
particularly  the  storm  decay  model,  provide  a  description  of  the  wind 
field  that  is  probably  less  accurate  at  sites  located  north  of  Cape 
Hatteras,  N.C.,  (mileposts  2200  or  more,   see  Figure  1).     The  numerical 
results  obtained  for  these  sites  should  therefore  be  viewed  very 
cautiously . 


24 


The  physical  model  of  the  hurricane  wind  speeds  may  be  refined  by  using 
a  detailed  representation  of  the  coastline  and  Its  adjacent  bodies  of 
water.     Such  refinements  may  result  in  an  Improvement  of  hurricane  wind 
speed  estimate  by  an  amount  of  the  order  of  5%  or  so. 

Additional  refinements  might  include  accounting  for  slower  storm  decay 
where  large  bodies  of  water  are  present  inland,  e.g.,  in  the  low- lying 
areas  of  Louisiana.     Additional  research  is  required  to  develop  Improved 
storm  decay  models  in  such  situations. 
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Facing  page:    Damage  caused  by  Hurricane  Camille. 
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5.  CONCLUSIONS 


In  this  paper,  estimates  are  presented  of  hurricane  wind  speeds  along 
the  Gulf  and  East  Coasts  of  the  United  States.     The  paper  describes  the 
sources  of  data,  the  probabilistic  models  for  the  climatological  charac- 
teristics of  hurricanes  and  the  physical  models  for  the  hurricane  wind 
speed  field  used  in  the  estimations.     Estimated  values  of  fastest-mile 
hurricane  wind  speeds  at  10  m  above  ground  in  open  terrain  at  the 
coastline  and  at  200  km  inland  are  given  for  various  mean  recurrence 
intervals  (Figures  1,  5,  and  6).     Fastest-mile  wind  speeds  at  distances 
inland  of  less  than  200  km  from  the  coastline  may  be  obtained  by 
linear  interpolation  between  these  values.     It  was  found  that  the 
estimated  hurricane  wind  speeds  are  best  fit  by  Weibull  distributions 
with  tail  length  parameters  y  >  A,  rather  than,   e.g.,  by  Extreme  Value 
Type  I  distributions.     Estimates  are  given  of  various  errors  inherent 
in  the  estimated  values  of  the  hurricane  wind  speeds.     The  confidence 
bands  for  the  estimates  were  found  to  be  of  the  order  of  at  least 
+  10%  at  the  68%  confidence  level.  It  was  also  noted  that,  owing  to 
the  possible  inapplicability  of  the  physical  models  used  in  this  work 
at  locations  north  of  Cape  Hatteras,  estimated  hurricane  wind  speeds 
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given  for  these  locations  should  be  viewed  with  considerable  caution. 
In  certain  special  cases  it  may  be  desirable  to  base  estimates  of 
hurricane  wind  speeds  on  more  elaborate  data,  physical  models,  and 
representations  of  the  coastline  than  those  used  in  this  work.  Never- 
theless, in  the  writers'  opinion,  the  results  presented  in  this  paper 
provide  a  rational  and  consistent  basis  for  making  realistic  decisions 
on  the  specification  of  design  wind  speeds  and  load  factors  in 
hurricane  -  prone  regions. 
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APPENDIX  A 


HURRICANE  WIND  SPEEDS  BLOWING  FROM  SPECIFIED  DIRECTIONS  FOR  lO-YR, 
50-YR,  and  lOO-YR  MEAN  RECURRENCE  INTERVALS 

Notes 

1.  For  key  to  mean  recurrence  interval  symbols,  see  bottom  of  next  page. 

2.  Ordinates  represent  estimated  fastest-mile  mean  speeds  in  mph 
(1  mph  =  0.4A7  m/s). 

3.  Abscissas  represent  clockwise  angles  of  directions  from  which  winds 
are  blowing  with  respect  to  north. 
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teristics of  building  elements  and  systems. 

Technical  Notes — Studies  or  reports  which  are  complete  in  them- 
selves but  restrictive  in  their  treatment  of  a  subject.  Analogous  to 
monographs  but  not  so  comprehensive  in  scope  or  definitive  in 
treatment  of  the  subject  area.  Often  serve  as  a  vehicle  for  final 
reports  of  work  performed  at  NBS  under  the  sponsorship  of  other 
government  agencies. 

Voluntary  Product  Standards — Developed  under  procedures 
published  by  the  Department  of  Commerce  in  Part  10,  Title  15,  of 
the  Code  of  Federal  Regulations.  The  standards  establish 
nationally  recognized  requirements  for  products,  and  provide  all 
concerned  interests  with  a  basis  for  common  understanding  of  the 
characteristics  of  the  products.  NBS  administers  this  program  as  a 
supplement  to  the  activities  of  the  private  sector  standardizing 
organizations. 

Consumer  Information  Series — Practical  information,  based  on 
NBS  research  and  experience,  covering  areas  of  interest  to  the  con- 
sumer. Easily  understandable  language  and  illustrations  provide 
useful  background  knowledge  for  shopping  in  today's  tech- 
nological marketplace. 

Order  the  above  NBS  puhlicaliviis  from:  Superintendent  oj  Docu- 
ments, Government  Printing  Office,  Washington,  DC  20402. 
Order  the  following  NBS  publications — FIPS  and  N BSIR's—from 
the  National  Technical  Information  Services,  Springfield,  VA  22161 . 

Federal  Information  Processing  Standards  Publications  (FIPS 
PUB) — Publications  in  this  series  collectively  constitute  the 
Federal  information  Processing  Standards  Register.  The  Register 
serves  as  the  official  source  of  information  in  the  Federal  Govern- 
ment regarding  standards  issued  by  NBS  pursuant  to  the  Federal 
Property  and  Administrative  Services  Act  of  1949  as  amended. 
Public  Law  89-306  (79  Stat.  1127),  and  as  implemented  by  Ex- 
ecutive Order  11717(38  FR  12315,  dated  May  II,  1973)  and  Part  6 
of  Title  15  CFR  (Code  of  Federal  Regulations). 

NBS  Interagency  Reports  (NBSIR) — A  special  series  of  interim  or 
final  reports  on  work  performed  by  NBS  for  outside  sponsors 
(both  government  and  aon-government).  In  general,  initial  dis- 
tribution is  handled  by  the  sponsor;  public  distribution  is  by  the 
National  Technical  information  Services,  Springfield,  VA  22161. 
in  paper  copy  or  microfiche  form. 


BIBLIOGRAPHIC  SUBSCRIPTION  SERVICES 


The  following  current-awareness  and  literature-survey  bibliographies 
are  issued  periodically  by  the  Bureau: 

Cryogenic  Data  Center  Current  Awareness  Service.  A  literature  sur- 
vey issued  biweekly.  Annual  subscription:  domestic  $25:  foreign 
$30. 

Liquefied  Natural  Gas.  A  literature  survey  issued  quarterly.  Annual 
subscription:  $20. 


Superconducting  Devices  and  Materials.  A  literature  survey  issued 
quarterly.  Annual  subscription:  $30.  Please  send  subscription  or- 
ders and  remittances  for  the  preceding  bibliographic  services  to  the 
National  Bureau  of  Standards,  Cryogenic  Data  Center  (736) 
Boulder,  CO  80303. 
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